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^ Abstract 

We examine in detail the diagrammatic mechanisms which provide the change of sign 
between the single transverse spin asymmetries measured in semi-inclusive deep inelastic 
scattering (SIDIS) and in the Drell-Yan process (DY). This asymmetry is known to arise 
^r^ due to the transverse spin dependence of the target proton combined with a T-odd complex 

(Nl phase. Using the discrete symmetry properties of transverse spinors, we show that the 

required complex phase originates in the denominators of rescattering diagrams and their 
respective cuts. For simplicity, we work in a model where the proton consists of a valence 
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^y-^ quark and a scalar diquark. We then show that the phases generated in SIDIS and in DY 

^-H originate from distinctly different cuts in the amplitudes, which at first appears to obscure 

^ the relationship between the single-spin asymmetries in the two processes. Nevertheless, 

further analysis demonstrates that the contributions of these cuts are identical in the 
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;_! leading-twist Bjorken kinematics considered, resulting in the standard sign-flip relation 
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between the Sivers functions in SIDIS and DY. Physically, this fundamental, but yet 
untested, prediction occurs because the Sivers effect in the Drell-Yan reaction is modified 
by the initial-state "lensing" interactions of the annihilating antiquark, in contrast to the 
final-state lensing which produces the Sivers effect in deep inelastic scattering. 



PACS numbers: 12.38.Bx, 13.88.+e 



1 Introduction 

1.1 Factorization and Final State Interactions 

The factorization picture of leading twist-perturbative QCD has played a guiding role in vir- 
tually all aspects of hadron physics phenomenology. In the case of inclusive reactions such 
as hadroproduction at large transverse momentum p + p — )■ H + X [Tll2], the parton model 
for perturbative quantum chromodynamics (pQCD) predicts that the cross section at leading 
order in the transverse momentum px can be computed by convoluting the perturbatively cal- 
culable hard subprocess quark and gluon cross section with the process-independent structure 
functions of the colliding hadrons and the final-state quark or gluon fragmentation functions. 
The resulting leading-twist cross section Enda/d^pHipp — ^ HX) scales as 1/p'^, modulo the 
DGLAP scaling violations derived from the logarithmic evolution of the structure functions 
and fragmentation distributions, as well as the running of the QCD coupling appearing in the 
hard-scattering subprocess matrix element. 

The effects of final-state interactions of the scattered quark in deep inelastic scattering are 
characterized by a Wilson line. Such effects have traditionally been assumed to either give 
an inconsequential phase factor or power-law suppressed corrections in hard pQCD reactions. 
However, this expectation is only true for sufficiently inclusive cross sections. For example, 
consider semi-inclusive deep inelastic lepton scattering (SIDIS) on a transversely polarized 
target i+p'^ — >■ H + i' +X. (For a review see [3|.) In this case the final-state gluonic interactions 
of the scattered quark lead to a pseudo-T-odd non-zero spin correlation of the lepton-quark 
scattering plane with the polarization of the target proton |1] which is not power-law suppressed 
with increasing virtuality of the photon Q^; i.e., it Bjorken-scales. This asymmetry is made 
experimentally explicit by studying the triple product of vectors Sp ■ {qx pn), where pn is the 
momentum of the hadron fragmented from the struck quark jet. Similar correlations can be 
found from target spin asymmetries, due to multi-photon exchanges [5j. 

A crucial fact is that the leading-twist "Sivers effect" j6l [7] is non-universal in the sense that 
pQCD predicts an opposite-sign correlation in Drell-Yan reactions relative to semi-inclusive 
deep inelastic scattering JHIIS]- This fundamental, but yet untested, prediction occurs because 
the Sivers effect in the Drell-Yan reaction is modified by the initial-state "lensing" interactions 
of the annihilating antiquark, in contrast to the final-state lensing which produces the Sivers 
effect in deep inelastic scattering. 

The calculation of the Sivers single-spin asymmetry in deep inelastic lepton scattering in 
QCD is illustrated schematically in Fig. [T| Although the Coulomb phase for a given partial 
wave is infinite, the interference of Coulomb phases arising from different partial waves leads to 
observable effects. The analysis requires two different orbital angular momentum components: 
5'-wave with the quark-spin parallel to the proton spin and P-wave for the quark with anti- 
parallel spin; the difference between the final-state "Coulomb" phases leads to a, S ■ {q x p) 
correlation of the proton's spin with the virtual photon-to-quark production plane |1|. Thus, 
as it is clear from its QED analog, the final-state gluonic interactions of the scattered quark 
lead to a pseudo-T-odd non-zero spin correlation of the lepton-quark scattering plane with the 
polarization of the target proton ^. The effect is pseudo-T odd due to the imaginary phase 
generated by the cut of the near-on-shell intermediate state. 

The S- and P-wave proton wavefunctions also appear in the calculation of the Pauli form 
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Figure 1: Origin of the Sivers single-spin asymmetry in semi-inclusive deep inelastic lepton 
scattering. 

factor quark-by-quark. Thus one can correlate the Sivers asymmetry for each struck quark 
with the anomalous magnetic moment of the proton carried by that quark [Til], leading to the 
prediction that the Sivers effect is larger for positive pions as seen by the HERMES experiment 
at DESY dl], the COMPASS experiment [El [131 [3 US] at CERN, and CLAS at Jefferson 
Laboratory [T6 | I17 | . 

The final-state interactions of the produced quark with its comoving spectators in SIDIS 
also produces a final-state T-odd polarization correlation - the "Collins effect", in which the 
Collins fragmentation function, instead of the Sivers distribution function, generates the spin 
asymmetry [T8| [T9] . This can be measured without beam polarization by measuring the cor- 
relation of the polarization of a hadron such as the A baryon with the quark-jet production 
plane [20l |2T| [22| [23] . Analogous spin effects occur in reactions in QED due to the rescattering 
via final-state Coulomb interactions. 

In principle, the physics of the "lensing dynamics" or Wilson-line physics [24, 25j underlying 
the Sivers effect involves nonperturbative quark-quark interactions at small momentum transfer, 
not the hard scale Q^ of the virtuality of the photon. These considerations have thus led 
to a reappraisal of the range of validity of the standard factorization ansatz. As noted by 



Collins and Qiu [2S], the traditional factorization formalism of perturbative QCD fails in detail 
for many hard inclusive reactions because of initial- and final-state lensing interactions. For 
example, if both the quark and antiquark in the Drell-Yan subprocess q + q ^ fi~^ + fi" interact 
with the spectators of the other hadron, then one predicts a cos 2(j) sin^ 9 planar correlation 
in unpolarized Drell-Yan reactions [27J. (Here (p is the angle between the niuon plane and 
the plane of the incident hadrons in the lepton pair center of mass frame, while 6 is the 
angle between the momentum of one of the muons and the line of flight of partons in the 
same frame [27> .28j.) This "double Boer-Mulders effect" [19j can account for the anomalously 
large cos 20 correlation and the corresponding violation [271 [29] of the Lam- Tung relation 
[301 [3T] for Drell-Yan processes observed by the NAlO collaboration |28|, and the azimuthal 
angle dependence of di-jet production in unpolarized hadron scattering (3^. Such effects again 
point to the importance of the corrections from initial and final-state interactions of the hard- 
scattering constituents, which are not included in the standard pQCD factorization formalism. 
The final-state interactions of the struck quark with the target spectators [33j also lead to 
diffractive events in deep inelastic scattering (DIS) at leading twist, such as £-|-p — )■ i' + p' + X, 
where the proton remains intact and isolated in rapidity; in fact, approximately 10 % of the deep 
inelastic lepton-proton scattering events observed at HERA are diffractive jSHES]- The presence 
of a rapidity gap between the target and diffractive system requires that the target remnant 
emerges in a color-singlet state; this is made possible in any gauge by the soft rescattering 
incorporated in the Wilson line or by augmented light-front wavef unctions. 

1.2 The Sign of the Sivers Effect 

As we have emphasized in the above, the sign reversal predicted by QCD between the Sivers 
functions jHl [71 [9] in semi-inclusive deep inelastic scattering (SIDIS) and in the Drell-Yan 
process (DY) [8, 9J is in contrast to the standard expectations of coUinear pQCD/parton model 
factorization. This reversal results in different signs predicted for the single transverse spin 
asymmetry (SSA) in SIDIS on a transversely polarized proton versus that in the DY process 
involving a transversely polarized hadron. The sign difference of the Sivers functions is a 
consequence of their being pseudo T-odd, with the SIDIS Sivers function related to the DY 
Sivers function by a time-reversal transformation |8l [361 E3|- The essential physics of the sign 
change is illustrated for QED in Fig. [2j The final-state interaction in SIDIS is attractive whereas 
the initial-state interaction is repulsive in DY. However, the detailed analysis involves many 
subtle aspects due to intermediate Glauber-type cuts related to unitarity. These considerations 
are worked out in detail in this article. 

In general, the generation of a SSA requires two ingredients: (i) the dependence of the cross 
section on the transverse polarization, and (ii) a relative complex phase between the amplitudes 
corresponding to different orbital angular momenta [351 [311 SO]- In the original argument given 
by Collins [8j the sign flip between SSA in SIDIS and DY arises due to a reversal of the direction 
of the Wilson line in the definition of the Sivers function under the application of T-reversal. It 
is interesting and important to clarify the relation between this time-reversal argument and the 
diagrammatic conditions (i) and (ii) stated above. While it is well-known that the propagation 
of a high energy quark or gluon can be approximated by a Wilson line along the corresponding 
light cone direction, the connection between the Wilson lines and Feynman diagrams can be 
non-trivial [41] . especially in the SSA case, where we are interested in the imaginary part of 
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Figure 2: The underlying physics of the sign change of the lensing interactions in SIDIS and 
DY reactions for QED. Dashed hues denote the photon propagators. The final-state interaction 
in SIDIS is attractive whereas the initial state interaction is repulsive in DY. 



the scattering amplitude due to the condition (ii) above. Establishing a clear connection in the 
Feynman diagram language between the spin asymmetries in SIDIS and DY will greatly aid 
future (model) diagrammatic calculations for the processes. 

Although the time-reversal symmetry argument allows one to confidently predict the sign flip 
between the Sivers functions in the two processes, a detailed diagrammatic model calculation 
demonstrating the origin of this sign reversal appears to be lacking in the literature. A model 
calculation for the SSA in SIDIS was first constructed in [9] by some of the authors; however, 
an analogous calculation for DY outlined in [4j assumed that the relative phase arises in DY 
due to putting the same propagators as in SIDIS on mass shell. As we will show below, one 
can think of the relative phase from the condition (ii) as a cut through an amplitude (or a 
complex conjugate amplitude), in addition to the standard final-state cut one obtains when the 
scattering amplitude is squared. In terms of the cuts, the calculation of [4j assumed that the 
cut generating SSA in SIDIS was equivalent to the cut in DY. While this assumption allows 
one to obtain the correct answer with the sign flip in the Sivers function between SIDIS and 
DY, as we shall discuss below, it is not easy to justify. In fact, the explicit calculation of the 
SSA in the SIDIS and DY processes given in this paper demonstrates that the cuts are, in 
fact, different, as illustrated in Figures |4] and [6j In Feynman diagram language this means that 
different propagators are put on the mass shell in order to extract the asymmetry in the two 
processes. In light-front perturbation theory (LFPTH) |^ H3] language this corresponds to 
putting different intermediate states (that is, intermediate states involving different particles) 
on the P~ LF energy shell. 



One may then wonder how the simple sign-flip relation between the Sivers functions in 
SIDIS and DY would arise from Feynman diagrams with different cuts (or LFPTH diagrams 
with different intermediate states put on energy shell). While further clarification of the dia- 
grammatic origin of the sign flip is left for future work, we shall show here by performing an 
explicit calculation that in the limit of high energy s and high photon virtuality Q^, the simple 
sign-flip relationship is preserved for SIDIS and DY Sivers functions, despite the different cuts 
in the diagrams in the two cases. We thus are able to construct a diagrammatic description of 
the sign reversal. 

The paper is structured as follows. In Section [2] we show how the transverse polarization 
dependence always comes into the amplitude squared with the imaginary factor of i. Since 
the cross section contribution has to be real, this implies that, for the spin asymmetry to be 
non-zero, another factor of i has to be generated elsewhere in the amplitude squared. We thus 
demonstrate the need for a phase difference between the amplitude and the complex conjugate 
amplitude contributing to the SSA [391140) . Our analysis provides an explicit verification of the 
analytical properties of the complete amplitudes. 

In Section [3] we perform a model calculation of the spin asymmetry in the SIDIS and DY 
cases using a Feynman diagram approach. We employ a model in which the proton is made 
out of a valence quark and a scalar color-charged diquark: this model was used previously in 
in [9]. Instead of calculating the Sivers functions directly, we consider the 'y* + p'^ ^ q + X 
process in place of the SIDIS Sivers function, and the g + p^ — )■ 7* + X process in place of 
the DY Sivers function. It can be shown that the relation between these processes and the 
corresponding Sivers functions is the same in both cases (see j27j for an example of how to 
establish the connection): hence, the sign reversal of the SSAs in the 7* + p^ — )■ q + X and 
g + p^ — )■ 7* + X processes is equivalent to the sign flips in the corresponding Sivers functions. 

Performing an explicit calculation in Section |3] we demonstrate that, in the framework of 
the model considered, the SSAs in SIDIS and DY arise from different cuts, as shown in Figs. |4] 
and |6j Nevertheless, the resulting spin asymmetries only differ by a minus sign, in agreement 
with the general arguments based on time-reversal anti-symmetry [H!^]- We demonstrate how 
the same conclusion can be achieved in a LFPTH calculation in Section HI while different LF 
energy denominators are put on energy shell in SIDIS and DY, the resulting SSA in DY is 
simply a negative of that in the SIDIS case. Section |5] summarizes our main results. 

2 Coupling of Transverse Spinors to the Complex Phase 

The discrete symmetries of the Dirac equation, C, P and T provide strong constraints on the 
form of the spinors and their matrix elements. In particular, the PT properties of transverse 
spinors bring in a unique coupling of the transverse spin dependence of a matrix element to 
its real and imaginary parts. This gives rise to the conclusion that the single transverse spin 
asymmetry requires the presence of a complex phase aside from the spinor matrix elements 
themselves. 

Below, we explicitly construct the transverse spinors and the identities they satisfy cor- 
responding to C/P/T transformations. We demonstrate that, as a consequence, when two 
spinor matrix elements are multiplied together (as in the numerator of a scattering amplitude 
squared), the part which depends on the transverse spin is pure imaginary. 



We then define the asymmetry observable A/^r and use the knowledge of the C/P/T con- 
straints to identify which types of diagrams can contribute to the asymmetry. We find that 
only interference terms can contribute to An, and furthermore only the imaginary part of the 
interference term (aside from the spinor products) contributes. We derive an expression for the 
spin-difference part of the amplitude squared in terms of the imaginary part of the rest of the 
interference term. 



2.1 Spinor Conventions and C/P/T Identities 

In this calculation we will work with the spinors defined in Ref. |l2]. These spinors correspond 
to the spin projection along the ±z axis of a particle with mass m and are expressed in light-cone 
coordinates p^ = p^^^ ±p^^\ For definiteness, working in the standard (Dirac) representation 
of the Clifford algebra, this spinor basis is: 
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Like any spinor basis for solutions of the Dirac equation, these spinors satisfy identities that 
embody the discrete C, P, and T symmetries of the theory. As can be explicitly verified from 
([I|, these spinors obey the identities 



C : -tl'VUp) = U^zip) 

PT : 7WVf/i.(p) = TU^zip) 

CPT : f/±,(p) = ±-f'V:^M, 



(2) 



where the final CPT identity combines the other two in a compact form. 

For our purposes, we are interested in transverse spin states. If we choose a frame in which 
the incoming polarized particle moves along the +z axis such that p± = 0±, then the spinors 
([I| become eigenstates of the helicity operator. We can form transverse spinors as done in [44j 
by taking linear combinations of ([I| to obtain the projection along, say, the x axis: 
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(3) 



where x = il is the spin eigenvalue along the x axis. Again, for a p± = 0± incoming particle, 
these spinors reflect a spin projection transverse to the beam axis; they are simultaneous eigen- 
states of the Dirac operator as well as the Pauli-Lubanski vector Wi, where W^ = —\efj,ypaS"'^p". 



Using ([2]) in (pi) gives the somewhat different C/P/T identities satisfied by the transverse 
spinors: 

C: -^7X*(P) = ^-x(p) (4) 

CPT : U^ip) = -xi'Vxip). 

Employing and generahzing Q allows us to write a complete set of identities for any transverse 
spinor matrix element: 

C: V^ik)j'^'---r"V^ip) = [V^iph''"---Y'Vyik)Y (5) 

C: U^ik)r'---r"V^ip) = [V^iph''"---Y'U^'{k)Y (6) 

CPT: V^ik)r'---r"U^{p) = Xx[V-^'ik)Y'---r-U^Ap)Y (7) 

u^ik)Y'---Y''U^ip) = xx[u^^'ik)Y'---Y''U^xip)Y- 

These identities allow us to explicitly determine rigid constraints on the form of any trans- 
verse spinor product. In particular, consider the parameterizations of both classes of spinor 
products: 

F^,(fc)7'^^ ■ --r-Uxip) ^ S^Aa{k,p) + xaik,p)] + 5^,^AKk.p) + xh\k,p)] (8) 

UAk)l^' ■ ■ ■ I'-U^ip) ^ SxAcik,p) + xc\k,p)] + S^,.Adik,p) + xd'{k,p)]- 

applying ([T]), one readily concludes that C/P/T constraints imply that: 

• a, b', c, and d' are real-valued. 

• a', b, c', and d are pure imaginary. 

Furthermore, this implies that if we multiply any two of these spinor matrix elements and 
sum over one of the spins (%'), e.g.. 

Ex' [Vx'ik)Y'---r"U^ip)] [Ux'ik)Y'---r"U^ip)r= (9) 

= [ac* + a'ic'y + bd* + b'id')*] + X [a{c'y + a'c* + b{d'y + b'd*], 

V V 

real imaginary 

we find that the product of two transverse matrix elements naturally partitions into a spin-even, 
real contribution, and a spin-odd, imaginary contribution. 

Thus in particular, the spin-dependent part of any product of two transverse matrix elements 
(say Si{x) and 5*2 (x)) is always pure imaginary: 

Si{x)s;{x) - s,{~x)s;{-x) = -[SKx)S2{x) - Sl{-x)S2{-x)]. (10) 



2.2 Relation of the Asymmetry to the Complex Phase 

The single transverse spin asymmetry An is an observable defined as 

_ t/a^(gl) - da^{q±) da^{q^) - da^{~q±) 
^ = W^JZ = W^JZ 'y^^) 

where da{q±) stands for the invariant cross section, e.g. ^^'^^ , for the production of a particular 
tagged particle with transverse momentum gl coming from scattering on a target with trans- 
verse spin 'I', I. Also daunp represents the normal unpolarized cross section ^{da'^ + da^). From 



(11), we see that the asymmetry An is proportional to the difference between the amplitude- 
squared l^p for X = +1 cind x = ~1- We denote this spin-difference amplitude squared as 
A\A\^: 

An (X \A\\x = +1) - \A\\x = -1) = A|^p. (12) 

Now let us identify the types of diagrams from which An can arise. Suppose there is a 
contribution from the square of an amplitude A{x) = F S{x) consisting of a spinor product 
S{x) which depends on the transverse spin eigenvalue x and a factor F coming from the rest of 
the diagram. Then the contribution of the square of the amplitude A to the asymmetry would 
be 

A\A\' ^ \An+l) - \An-l) (13) 

= \F'\ [\S{+1)\'-\S{-1)\']. 



But the constraints on the spinors due to C/P/T (10) imply, for Si = S2 = S, that 

[\S{+1)\'-\S{-1)\']=0. (14) 

This is easy to understand mathematically: the spin-dependent part must be pure imaginary, 
but any amplitude-squared is explicitly real. Hence, any amplitude squared is independent of 
X and cannot generate the asymmetry; An can only be generated by the quantum interference 
between two different diagrams. 

Thus at lowest order in perturbation theory, the asymmetry could be generated by the 
overlap between an 0{as) one-loop virtual correction and the Born-level amplitude. So let us 
consider a similar exercise to determine the contribution to A|^p from this 0{as) correction 
(as compared to the Born- level amplitude squared). Let us write the tree- level amplitude ^(0) 
and the one-loop amplitude ^(1) as 

Ao)(x) = F.s^ix) 

where the factor Si includes all momentum and spin-dependent numerators, and the factor Di 
contains all the propagator denominators. At 0{as) (as compared to the Born-level amplitude 
squared) the spin-difference contribution is 

A\A\'' = Ai)(+iMm(+i) + ^;i,(+iM(o)(+i) - (x^-x) (16) 



But from the C/P/T constraints (10), we see that the numerator of (16) is pure imaginary; 
thus 
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Thus we conclude that the spin-dependent part which contributes to the asymmetry comes 
only from the imaginary part of the remainder of the interference term, aside from the spinor 
matrix elements. This is easy to understand mathematically: if the spin-dependent part of 
the spinor matrix elements is pure imaginary, then it must multiply the imaginary part of the 
remainder of the interference term to generate a real contribution to the asymmetry. This 
imaginary part of the rest of the amplitude interference term gives the complex phase that is 
required by C/P/T to generate the asymmetry A^; it is not simply the imaginary part of any 
one diagram, but rather a relative phase between the tree-level and one-loop amplitudes. If 
there is no relative phase present in the pre-factors, e.g. Im(FiF2 ) = 0, then the imaginary part 
comes from the denominator of the loop integral Di{k). In that case, taking the imaginary part 
corresponds to putting an intermediate virtual state on energy-shell in light-front perturbation 
theory. The imaginary part generated this way was discussed in |1] and [9] as a source of the 
Sivers-type asymmetry. 

3 Model Calculations with Feynman Diagrams 

Following [^ and [U], we employ a toy model of a point-like proton. To model the parton 
distribution function, we introduce a Yukawa-type coupling between the proton field, the quark 
field, and a scalar "diquark" field. This corresponds to an interaction term in the Lagrangian of 
the form +G{ip ipq + iljq'ipp)4'qq- The scalar diquark field couples to gluons by the rules of scalar 
QED, using the same covariant derivative iD^j, = id^ + gA°jT'^ as the fermions. We use g to 
represent the QCD coupling of the quark and diquark and e/ to represent the electromagnetic 
charge of the quark of flavor /. 

In the next sections we will explicitly calculate the spin-difference amplitudes at 0{as) 



defined in (16) for deep inelastic scattering and for the Drell-Yan process. We will demonstrate 



in detail the emergence of the predicted minus sign difference between the Sivers asymmetries 
in the two processes. Throughout the calculation, we will work with the light-cone coordinates 
and corresponding metric 
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Additionally, we will drop the quark mass m wherever it occurs, but we will keep the proton 
mass M and the mass A of the scalar diquark field. When it is necessary to approximate the 
kinematics, we will work in the limit 
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Figure 3: Diagrams for the 'y* + p^ ^ q + X SIDIS amplitude at one-loop order (A) and tree- 
level (B). The incoming solid line denotes the transversely polarized proton, which splits into 
a quark (outgoing solid line) and a diquark (dashed line). 



3.1 Semi- Inclusive Deep Inelastic Scattering 

For deep inelastic scattering, we consider the scattering of a virtual photon with virtuality 
g^ = —Q^ on a transversely polarized proton with transverse spin eigenvalue x- At lowest 
order, this process produces a quark and diquark, as shown in Fig. [3j 

To begin, let us establish the kinematics. Following [1], we work in the Drell-Yan-West 
frame which is coUinear to the proton {pj_ = 0^) and boosted such that g+ = exactly. In 
this frame, then, the photon's virtuality comes from its transverse components: Q"^ = q]_. We 
define the longitudinal momentum fraction exchanged in the t-channel as A = r^ jp^ . Then 
momentum conservation and the on-shell conditions for the proton, quark, and diquark fix r~ 
and q~: 
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With these kinematics, we can write down the one-loop amphtude shown in Fig. pi (A) as 
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where the longitudinal momentum fraction in the loop is k~^ = xp^ and Cp = {N^ — '^)/'iNc 
is the Casimir operator in the fundamental representation. Similarly, the tree-level amplitude 
shown in Fig. Is] (B) is 



e/G, 



A^''=^-^UAq + r)iJU^{p). 



(23) 



The lowest-order contribution to the spin-difference amplitude squared comes from the over- 



lap of these diagrams and, in particular, the imaginary part of the denominators (cf. Eq. (17)): 



AU 



DIS\ 



2i 



g^ejG^Cp 
2(27r)V(j9+)3 



dxd'^k 



x^(x — A)(l — x) 
1 



Im < dk 



X 




(24) 



k- + q- - 



{q±+k±)^-ie 



xp'^ 



h- — r- ^ {k±-r±y-ie 
'^ ^ (x-A)p+ 



X J2[UxiP)frxUx'iQ + r)U^iq + 
where we sum over the spin of the outgoing quark and use 



_^" -P'+ (l-x-)p+ 



(25) 



We will provide a justification for summing over photon polarizations in the two paragraphs 



following Eq. (|62|). Performing these sums and simplifying the expression gives 

dxd'^k 



A\A 



2ig^e]G^CF 



DIS\ 



X 



{2^Yr^{p+f J x2(x-A)(l-a;) 
f/x(p)/(^ + mf-t- m + f)tUx{p) -ix^ -X) 



(26) 
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where the imaginary part that is essential for generating the asymmetry comes from the ex- 
pression 



X = Im< 



i dk 



k~ - 



xp~' 






k — r 



(fej 



J^ 



{x-A)p+ 



k — p + 



k'j_+\'2-ie 
{l-x)p+ 



(27) 



Notice that the numerator in (26) containing the Dirac matrix element is ^"-dependent, 
such that the dk~ integration contained in X from Eq. (27) applies to it too. Superficially 
the numerator of (26) could scale as {k~Y at large k~ , which would endanger convergence; 
however, it actually only scales as {k~y since (7"*")^ = 0. Thus the k~ integral scales at most 
as dk~ /{k~y, which converges and allows us to close the contour in either the upper or the 
lower half-plane. 



In addition, we will demonstrate below that in the kinematic limit at hand given by Eq. (19) 



the leading contribution to the Dirac matrix element in Eq. (26) is, in fact, k -independent 



We will, therefore, proceed under the assumption that this is the case and that all the k 



dependence in (26) is contained in the integrand of Eq. (27), evaluating the integration in (27) 
separately. 



The imaginary part in (27) comes from the denominators, which corresponds to putting two 



of the loop propagators on shell: one occurs from the residue of the dk~ integral and the other 
occurs by taking the imaginary part. However, there are strong kinematic constraints that 
restrict which combinations of propagators can go on-shell simultaneously. We are working in 
the limit of massless quarks, and 1 f-)- 2 processes for on-shell massless particles are forbidden by 
four-momentum conservation; cuts corresponding to such processes will explicitly be impossible 
to put on shell. Other cuts correspond to spontaneous proton decay; proton stability against 
decay through various channels must be imposed by hand, resulting in kinematic constraints 
on the masses of the proton and the scalar. 

Eq. (27) is evaluated in Appendix lAJ The result reads 



X 



27r2AV 




6 


( 

X — 
\ 


i + 2^^^-(^^r''^)) a] 




Q' 


p- 


kl+x^ 

{l-x)p+ 


k] 1 

xp+ 




kl+X2 
P ^ (l-a:)p^ 


(fcx-r±)2] 
- "•" {A-x)p+ 



(28) 



Substituting this expression back into (26) and integrating over the delta function which sets 
a; ~ A gives 



A\A 1^- ,«Vr^r iAi^^^\ f d'k 

A \Anis\ - -^g ep C, (g2(,-.+«2)J / (^ 



(29) 



f/. 



X 



+ 



-^ 



+ f)Wx{p)-ix^-x) 



(fc± - r^Yikl + a2) 



where the mass parameter 0? that regulates the infrared divergence 

a^ = A (A^ - (1 - A)M2) > 



(30) 
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is ensured to be positive definite by the proton stability conditions (A3) and (A7), and we have 
used 

Note that making this cut has fixed the loop nionientuni k'^ to be 

M2 A^2 ^ ^2 



kf" 



Ap^ 



p-\ 



[1 - A)p+ ' 



A;, 



(32) 



Next we need to evaluate the numerator of (29) by computing the difference between the 
matrix elements: 



ndis = u^{p)f{f +mf-t- m + f)tUxip) - (X ^ -X)- 

The momenta involved in this spinor product obey the scale hierarchy 

0(Q) Oil.) 0{^VQ) 



(33) 



p~^,r'^,k^,q ,\q±\ > \r±\,\k±\,M, \ ':$> p ,r ,k , (34) 

with the dominant power-counting of the spin-dependent part of the Dirac matrix element 
being 0{Q^1-'^). Evaluation of N^is in Eq. (33) in the kinematics of Eq. (34) is somewhat 



involved: after some algebra one can show that there are three classes of Dirac structures that 
give a contribution of the leading order, (^(Q^-L^); all three involve taking the 0{Q) momenta 
from the middle three gamma matrices: 



N, 



DIS 



_ [(2p+ -k+- r+){q-y] U^ip) f^+^-^+^ U^{p) - ix 
[(1 - A)p+(g-)2] U^{p)f^^^U^{p) - (x ^ -X). 



-x) 



(35) 



The three variations consist of taking 7" for both f and |^, taking 7" for one and 7_l for the 
other, or taking 7_l for both. 

In the first case, if we take 7" for both f and |^, we obtain 



U^{p)fl^l/^U^{p) ^ -A2(p+)2f/^(p) 7-7+7- f/x(p) 

= A\p+fu^{p)ru^{p), 



(36) 



but U^{p)'~) U^{p) = —f- is spin- independent and cannot generate the asymmetry. Similarly, 
if we take 7_l for both f and ft, we obtain 

U^{p)tl^l^U^{p) ^ r\k{U^{p)Y±l'-liU^ip), (37) 

but U^{p) 7^7"*" 7;^ U~^{p) = 2p+(5*-' is also spin-independent and cannot generate the asymmetry. 
However, if we take one each of 7_|_ and 7", we obtain 



f/^(p)/7+|^f/^(p) ^ --Ap 






(38) 
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We can further simplify this expression by using the Dirac equation 



Q = {f-M)U^{p) 



1 



rP 7 + 



M^ 



-r 



-M 



uAp) 



(39) 



2' ' 2p^ 

to rewrite the action of 7" in terms of 7+ and M. Since (7"'')^ = 0, this simphfies (38) to 

UAp)tl^^U^ip) ^ -MA{ki-ri)U^{p)^^Y±U^ip), (40) 

and U^{p)'y~^'y]_ U^{p) = 2ixp^S^'^, which is spin-dependent and generates the asymmetry. Al- 
together this gives 

Nois = -2zxA(l-A)(p+)2(g-)2M(fcf -rf)-(x^-x) 



(41) 



-4z 



A 



(2) J2h 



Q^'M^k'l' -rf) 



so that the spin-difference matrix element is pure imaginary, as was proved in (10) 
Substituting this result back into (29) gives (cf. Eq. (31) in [27|j) 



A|^ 



DIS\ 



-Ag'eJG'Cp 
( 

IT 



AfQ^M 



f^: + a'^ 



(Pk 



,(2) _ ^(2) 



(2vr)2 ik, 



TL? (kl + a2) 



A)^ 



Q^Mr 



(2) 



f^iif^i + a? 



In 



r^^ + a? 



where the d?k integral is performed using Feynman parameters obtaining 



(Pk 



k 



(2) 



,(2) 



(2) 



(27r)2 {k^ _ f^)2 (^2 + ^2) 



1 r___ 



In 



r^^ + 0? 



(42) 



(43) 



Eq. (42) is the final expression for the spin-difference amplitude squared for deep inelastic 



scattering. We would next like to compare this expression with the result for the Drell-Yan 
process. 

Before we do that, let us stress once again that the asymmetry in the SIDIS case arises from 
the contribution of the diagram in Fig. [3] (A) with the (g + k)- and [p — fc)-lines (corresponding 
to the lines labeled ® and (D in Appendix lAl) which are put on mass-shell. It is this and only 
this contribution that gives the imaginary phase needed for the asymmetry in SIDIS. This fact 
becomes more apparent if we diagrammatically represent putting the (g + k)- and {p — A;)-lines 
on mass-shell by a cut, as shown in Fig. |4j In Fig. |4]we show the amplitude squared which we 
have just calculated, with the longer cut representing the true final state of the process, and 
the shorter vertical cut line representing the imaginary phase generating the asymmetry. The 



shorter cut follows the standard Cutkosky rules [46j, with the caveat stressed above in Sec. 2.2 
that it should not be applied to the spinor matrix element; that is, the shorter cut applies to 
the denominators of the propagators only, as if we are evaluating the diagram in a scalar field 
theory. Using the Cutkosky rules one can clearly see that this is the only way the shorter cut 
line can be placed in the diagram, as all other cuts would lead to various prohibited 1 — )■ 2 or 
2—7-1 processes, including proton decay. Thus Fig. |4] demonstrates that the imaginary phase 
needed for the single-spin asymmetry arises only in diagrams where it is possible to place a 
second cut. We will make use of this result in the analysis of the Drell-Yan process below. 



^As noted in Ref. [45 , there should be an additional overall minus sign in front of Eq. (21) of Ref. |4], also 
in front of Eq. (31) of Ref. [9_ and Eqs. (31,33,36) of Ref . [27]. 
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Figure 4: Diagrammatic representation of the origin of complex phase leading to the single-spin 
asymmetry in SIDIS. The longer cut denotes the final state of the process, while the shorter 
cut demonstrates the origin of the phase needed for the asymmetry. 

3.2 Drell-Yan Process 

We now perform a similar calculation for the Drell-Yan process in the same model consid- 
ered above for deep inelastic scattering. We will consider the scattering of an antiquark on a 
transversely-polarized proton with transverse spin eigenvalue x that produces a virtual photon, 
which then decays into a dilepton pair with invariant mass (^ = Q"^ ■ This process is shown in 
Fig. [5] at the level of virtual photon production: g -|- p^ — )■ 7* -|- X. 



(^) 



X 



X 



q — r 



q — k 




WWW 



p 



p — k 



p — r 



(B) 



X 



X 



rvwwv ^ 



p 



p — r 



Figure 5: Diagrams for the q + p'^ —>'-/* + X DY amplitude at one- loop order (A) and tree- 
level (B). The incoming proton and anti-quark are denoted by the lower and upper solid lines 
correspondingly, with the outgoing diquark denoted by the dashed line. 

Following [9j, we work in a generic frame coUinear to the proton (pj_ = 0_l). We define the 
longitudinal momentum fraction of the photon to be /3 = q^ /p^ and the momentum fraction 
exchanged in the t-channel to be A = r^/p'^. As before, four-momentum conservation and the 
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on-shell conditions fix r and q to be 

p+ (1 - A)p+ 



P 



{p — r)~ 



(44) 



(g — r) + r" 



(gl - r^f 
(/3 - A)p+ 



+ r" 



^ g^ - 2gl • rl ^ ^ 



(/3 - A)pH 



jji 



In this frame, tlie virtual plioton's large invariant mass Q^ comes in part from its transverse 
components and in part from its longitudinal components: 



Q^^q' = ^p+q- ^ 



Qi 



A 



This allows us to approximate q as 



Q 



AJ9^ 



+ 



/3-A 
Ql 



qi + O 



Q 



pi 



(45) 



(46) 



which agrees with the corresponding expression (20) for DIS to leading order in Q^. The 
kinematics can be summarized as 



M2 -. 

M / « + iq±-r±? ^ M' 
1^ = {t^P , —o — W^ZZ + -^ 



(47) 



(/3 - A)pH 



rl + A^ 
(1 - A)p^ 



g± 



+ M^ rj + A^ ^ 
^ ' p+ (1 - A)p+ ' ""^ 

Notice that, in this frame, the on-shell conditions for the antiquark, scalar, and dilepton pair 
imply that (/3 — A) > 0, /3 > 0, and (1 — A) > 0. Additionally, to leading order, the positivity 



constraint on q (46) implies that A > 0, and we can choose our frame such that /3 < 1, 



although this is not strictly necessary. Altogether, this gives a hierarchy of the fixed scales to 
be < A < /3 < 1. 

With these kinematics, we can evaluate the one-loop amplitude shown in Fig. [s] (A) as 



A{ 



j^Y _ ig^efGCp 



(Tk, 



VAq-r){2f-]^-m-(^)i\mx{p) 



(27r)4 J [k-i + ie][{k - qf + ie][{k - ry + ie][{p - k)"^ - X^ + le] 

-ig'^CfGCp f dx dk^ d'^k 



(48) 



X 



2(27r)4(p+)3 7 x{x - /3){x - A){1 - x) 

VAq-r){2f-t-m-<l)i\mx{p) 



k- 



k^—ie 
xp+ 



k — q' 



(x'-/3)p+ 



U- _ r- — ik±-f±)^-ie 
'^ ' {x-A)p+ 



k — p 



k^+X^-ie 
{l-x)p+ 



where x = k~^ /p~^ is the longitudinal momentum fraction in the loop. Similarly, the tree-level 
amplitude shown in Fig. [s] (B) is 



•^0 



'-i^VAq-rWUAp)- 



(49) 
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This allows us to calculate the spin-difference amplitude squared following (17) as 



AU, '2 



^DY\ 



2i 



_2(27r)4(p+)3r2 



dxd k 
x(x-/3)(x- A)(l -x) 



Im 



dk 



k- - 



fc2. 



xp~' 



X 















1 








k- 


-q- 


{k±-q±y-ie 
{x-P)p+ 




k- 


— r^ 


(fc^-rj^)2-ie 
(x-A)p+ 




k- 


P' + (1-.)P+ J 



(50) 



X',A 



— r 



*JuM-{x^-x 



where we sum over the spin of the incoming antiquark and use Eq. (25). (See the discussion 



following Eq. (62) for a justification of summing over photon polarizations.) Performing these 

dxd'^k 



sums and simplifying the result gives 

2ig^e)G^CF 



/\\A 



DY\ 



X 



(27r)V(p+)3 7 x(a;-/3)(x-A)(l-x) 
X \U^{p) f{$ -(f){2f-l^- f){(^ - f)t U^ip) - (x ^ -x)] 

where the imaginary part necessary for the asymmetry is generated by 



(51) 



Im< 



k- - 



xpf 



k — q 



(x-P)p+ 



i dk 
k~ — r' 



{x-/S.)p+ 



k — p 



kl+X'^-ie 
{l-x)p+ 



(52) 



As before, the imaginary part (52) corresponds to putting two of the loop propagators on-shell 



simultaneously: one from performing the k integral and another from taking the imaginary 
part. The propagators that can be simultaneously put on-shell are strongly constrained by 



the kinematics and by the requirement of proton stability. The expression (52) is evaluated in 
Appendix |B] yielding 



I 



27r2A(/5 - A)p+ 




5 


X- (l + 2 


Q' ) \ 




Q' 


[xp+ 


— r~ 


1 (fex-rx)2" 

^ {A-x)p+\ 




[xp+ P' + 


kl+\2 " 
{l-x)p+\ 



(53) 



which we can substitute back into (51 ). Integrating over the delta function sets x ~ A, giving 

d'^k 



A\Ady\^ = -tg^eJG^Cp 



A(l-A) 



X 



Q^{ri + a'^)J J (27r)2 
UxJP) fit -i){2f-1,- f){^ - f)i U^jp) -ix^ -X) 
ik± - r-l)2 {kl + a2) 



(54) 



where we have again employed (30) and (31 ), since we have established that the proton stability 



constraint (|A3|) is still valid for the Drell-Yan process. In performing the longitudinal integrals, 



we have fixed the loop momentum k^ to be 

kf" = I Ap 



Ap^ 



(55) 
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Comparison of (29) with (54) shows that the only difference between the two processes 



occurs in the numerators, rather than in the denominators. The essential difference in the 
numerators is the reversal of the intermediate (anti)quark propagator from (^+$ in deep inelastic 
scattering to |^ — ^ in the Drell-Yan process. We will return to this point later in the analysis 
of the results. 

Next we need to evaluate the spin-difference matrix element appearing in the numerator of 



(54): 



N, 



DY 



uAp) 



f)tUx{p)-{x-^-x)- 



(56) 



The momenta obey the same scale hierarchy (34) as in deep inelastic scattering, with the 



addition of q^ as a scale at 0{Q) in our frame for Drell-Yan. The other momenta can differ 
from their values in DIS by factors of 0{1), but the power-counting is the same. Again, the 
dominant power-counting of the matrix element is (9(Q^±^), which only arises from taking 



(^-^)(2|i-^- 



= -(l-A)(p+)(g-)V 



(57) 



so that 



Ndy = - [(1 - A)p+(g-)^] U^{p) i^H UAp) -ix^ -x)- 



Comparing (58) with (35), we see that 



Nnv = -N, 



DIS 



(58) 



(59) 



to leading order in Q, so we can immediately write the numerator for Drell-Yan using (41) as 

l-A^ 



Nny = +4i 



A 



Q*M{k][> - r 



(2) ^(2)^ 



(60) 



Substituting this back into (54) yields the same transverse momentum integral as in DIS, which 



we can evaluate using Feynman parameters to obtain the final answer 



A\A 



g'ejG'^CF 



DY\ 



TT 



1-A)^ 



Q^Mr 



(2) 



r^Af^i + a^] 



In 



f^^ + a^ 



(61) 



-A|^ 



DIS\ 



Thus we conclude that the spin-difference amplitude squared from the Drell-Yan process is 



exactly the negative of the that from deep inelastic scattering, (42). 



To obtain the single-spin asymmetry Aj^i one needs to divide A|^p for DY and SIDIS by 
twice the unpolarized amplitude squared (averaged over the incoming proton polarizations), as 



follows from Eq. (11). Both in the SIDIS and DY cases the unpolarized amplitude squared 



is dominated by the Born-level processes, with the amplitudes given in Eqs. (23) and (49) 



correspondingly. One can easily show that the squares of those amplitudes, averaged over the 



proton polarizations, are, in fact, equal, such that Eq. (61) leads to [81 [9] 



aDY 



-A 



DIS 

N ■ 



(62) 



The sign-reversal has also been derived in [S] for the SIDIS and DY Sivers functions [HI 
El El inS]- At the same time in our analysis we have studied the full SIDIS and DY processes 
(cf. [H |9]) instead of the corresponding Sivers functions. However, the conclusion (62) was 



reached using Eq. (25) above: in particular, note that we have reduced the SIDIS process to 
7* +p^ — > g + X scattering and have summed over polarizations of the incoming virtual photon. 
This is not an exact representation of the physical SIDIS process, since we have to convolute the 
hadronic interaction part of the diagram with a lepton tensor coming from the electron-photon 
interactions. Likewise, for the Drell-Yan process we have replaced the second hadron by an 
anti-quark, reducing it to the q + p^ ^ '-j* + X scattering. (Replacing the /"""/"-pair by a time- 
like photon is also an approximation, true up to an overall multiplicative factor due to current 
conservation.) The resolution for these questions is in the fact that, in the eikonal kinematics 



(19) considered, the single-spin asymmetries in the 7* + p''^ — )■ g + X and q + p^ ^ '-)* + X 



processes are proportional to the Sivers functions for SIDIS and DY correspondingly, as can be 
shown along the lines of [27J. 

Consider the quark correlator in a proton ^7\ |27] 



$,,(A,r-l;P,^) 



dx (fx± i(i 
(2 7r)3 



xp^ 



{P,S\tPj{0)Uij,{x+ = 0,x-,x^)\P,S), (63) 



where the quark has transverse momentum f± and the longitudinal momentum fraction A, 
the proton spin four- vector is S"^, while U is the gauge link necessary to make the object 
gauge- invariant. The correlation function $jj can be decomposed as [T9l [27] 



<^,,{A,rr,P,S) 



M 
2p^ 



1 



/i(A,r-l){- + — AVA,r-l)e^,,.7>^r^5?--gi.(A,f^)^7 



M M2 



M 



1 1 ^f^ pi' 

- j^ hiriA, f±) i a^^ 7^ ^^P'^ - jp ^ti^^ r±) i a^^ 7^ r^ p" + h^{A, f±) a^^ -^ 



(64) 



«j 



The Sivers function /£'"j,(A,rj_) can be singled out by extracting the spin-dependent part of 
$i,(7+),i, that is [27] 



^iAl 



)ji 



spin dependent 



M 



s'^ Sir\ 



fiA^,r^) 



(65) 



Comparing Eq. (65) to Eqs. (35) and (58) above (and comparing the latter two equations to 
Eq. (29) in |27j ) we see that both calculation performed here for the SIDIS and DY processes 
single out the corresponding Sivers functions. In fact our calculation is consistent with that 



performed in |27) . as can be seen by comparing Eqs. (42) and (61) to Eq. (31) in ^J. We can 
understand this consistency from the fact that the summing over the polarization of the virtual 
photon after Eq. (24) and (50) (see Eq. (25)) corresponds to replacing 7"*" in the left hand side 



of Eq. (65) by 7^7^7^ which is just —27+. Thus the contributions of both Eqs. (42) and (61) 



are proportional to the SIDIS and DY Sivers functions, calculated in the particular model for 
the proton considered here, with, as one can show, identical proportionality coefficients. (This 



point is strengthened further by noticing that the relation in Eq. (61 ) is only valid if one writes 
the spin-difference amplitudes in terms of Q^ and A = xp, as is proper for the distribution 
function like a Sivers function.) Therefore, the sign reversal in Eq. (62) is just an explicit 



manifestation of the sign reversal between the SIDIS and DY Sivers functions. 



19 




yVWWWAA/ 



Figure 6: Diagraniniatic representation of the origin of complex phase leading to the single-spin 
asymmetry in the Drell-Yan process. The longer cut denotes the final state of the process, while 
the shorter cut demonstrates the origin of the phase needed for the asymmetry. 



It is interesting to investigate the diagrammatic origin of the sign-flip in Eqs. (61 ) and (62). 



To do that we consider the diagram contributing to the single-spin asymmetry in the Drell-Yan 
process shown in Fig. |6j As follows from the calculation in Appendix |B| the asymmetry in 
the Drell-Yan case arises due to putting the (g — k)- and k-lines in Fig. [s] (A) (corresponding 
to lines ® and (2) in Figs. 13 and [l4|) on mass-shell: this is illustrated in Fig. |6| by the second 



rX—L. 



(shorter) cut, in analogy to Fig. |4j Comparing Figures |6] and 111 we see that the minus sign in 
Eqs. (61 ) and (62) arises due to the replacement of the outgoing eikonal quark in Fig. lilby the 



incoming eikonal anti-quark in Fig. [6j this is in complete analogy with the original Wilson-line 
time-reversal argument of Collins [8j (see also |36|). 

However, a closer inspection of Figures |4] and [6] reveals that the cuts generating the complex 
phase appear to be different: in Fig.lllthe (shorter) cut crosses the struck quark and the diquark 
lines, while in Fig. p\ the (shorter) cut crosses the anti-quark line and the line of the quark in 
the proton wave function. While we have already identified the outgoing quark/incoming anti- 
quark duality in SIDIS vs. DY as generating the sign flip, the fact that in the proton's wave 
function the diquark is put on mass shell in SIDIS and the quark is put on mass shell in DY 



makes one wonder why the absolute magnitudes of the asymmetries in Eq. (62 ) are equal. After 



all, different cuts may lead to different contributions to the magnitudes of the asymmetry. 



Ultimately the origin of Eq. (62) is in the fact that spin-asymmetry is a pseudo T-odd 



quantity and the Wilson lines describing the outgoing quark in SIDIS and the incoming anti- 
quark in DY are related by a time-reversal transformation [8j. However, in the diagrams at 
hand the origin of the equivalence of the shorter cuts in Figs. |4] and [6] is as follows. Consider the 
splitting of a polarized proton into a quark and a diquark as shown in Fig. [7j this subprocess 
is common to both diagrams in Figs. |4] and [6j The essential difference between Figs. |4] and |6] 
that we are analyzing is in the fact that in Fig. |4]the diquark is on mass shell, while in Fig. |6] 
the quark is on mass shell. 

Concentrating on the denominators of the quark and diquark propagators in Fig. [7] we shall 
write for the SIDIS case of Fig. |4] (quark is off mass shell, diquark is on mass shell) 






5\ k- 



M2 



fcj + A^ 

;i-A)p-t 



p+ {kl + a2) 



5{k-), (66) 



where we have used Eqs. (21), (34), and (30) along with x ~ A, and, in the last step, neglected 
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p-k 

Figure 7: Polarized proton splitting into a quark and a diquark, which is a part of the diagrams 
in both Figs. |4]and[6] 



all (9(X^ IQ) terms inside the delta-function since the numerator of the diagram does not 
depend in the exact value of k~ as long as it is small. 

A similar calculation for the Drell-Yan process from Fig. [6] (quark is on mass shell, diquark 



is off mass shell in Fig. 7| employing Eqs. (47) and (34) leads to 



(p - kf - X' 



5{e) 



p+ {kl + a2) 



6 k-- 



-1 



Ap-i 



p+{k] 



6{k- 



(67) 



We see that although the two contributions in Eqs. (66) and (67) are, in general, different, in 



the kinematics (34) they are apparently equivalent, leading to two different cuts in Figs. 4 and 



[6] giving the same-magnitude asymmetries. 

To complete this Section, let us note that, in the framework of the model at hand, there is an- 
other diagram in the Drell-Yan process which at first glance contains both the spin-dependence 
and a complex phase needed to generate the single-spin asymmetry. The diagram is shown in 
Fig. [8] with its contribution to the single-spin asymmetry denoted by the double-cut notation 
of Figs. |4]and|6} The potential contribution to the asymmetry arises due to a phase generated 
by the correction to the quark-photon vertex in Fig. |8j Note that an analogous graph cannot 
give an imaginary part in the case of SIDIS, since there the virtual photon is space-like. 




7AAAAAAA/AA4 



p-r 



q-r 



Figure 8: The potential contribution to the asymmetry in DY coming from the quark-photon 
vertex correction. 

We will also demonstrate that the contribution of the diagram in Fig. |8] to the single-spin 
asymmetry is zero. To do this one needs to evaluate the numerator of this diagram (minus the 
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spin-flip term): 

A 

= -8k-{q-r) U^{p) f{^ - ^)f U^{p) - (x ^ -X) 

= _8 A; . (g - r) [2r • (g - k)U^ip) fU^ip) - r'U^ip) (^ - ^) U^ip)] - (x ^ -x) = 0- (68) 

The zero answer results from the fact that, as can be checked explicitly, forward Dirac matrix 
elements of transverse spinors with a single gamma- matrix, i.e. expressions like U^{p) j'^ U^{p), 
are x-independent. Hence, the diagram in Fig. M does not contribute to the asymmetry. (In 



fact, the second line of Eq. (68) is proportional to the square of the Born term from Fig. M (B): 
as we show in Sec. 2.2| the square of the Born diagram cannot lead to a non-zero single-spin 
asymmetry. ) 

Finally, let us point out that in the calculation of the asymmetries in both SIDIS and DY, we 
have neglected diagrams in which the virtual photon couples to either the proton or the scalar 
diquark instead of the (anti)quark. These diagrams are necessary to ensure gauge invariance, 
but they are suppressed by powers of -L/Q, which allowed us to neglect them. 

4 Model Calculations with LFPTH 

In this section we show that we can re-derive the results obtained in Section [3] by using light- 
front perturbation theory (LFPTH) ^21 H3] . The LFPTH approach and covariant Feynman 
diagrams calculations are equivalent; however, we find it instructive to show the equivalence 
explicitly. 

4.1 SIDIS 

The LFPTH diagrams contributing to the single-spin asymmetry in SIDIS are shown in Fig. [9j 
Let us point out from the outset that the diagrams containing instantaneous terms do not 
contribute to the asymmetry in SIDIS and, therefore, are not shown in Fig. |9} This is clear 
from the calculation of the numerator in Sec. 13. 1[ The dominant contribution to the numerator 
comes from 7" in the quark-gluon vertex, which eliminated the possibility of the instantaneous 



gluon exchange contributing. The factor of 7"^ arising in the second line of Eq. ( 35 ) due to the 
photon-quark interactions eliminates the possibility of the instantaneous quark line exchanges. 

The numerators of the amplitudes calculated in the Feynman diagram approach of Sec. [3| 
consisting of the Dirac matrix element and vertex factors, are clearly identical to those one 
would find using LFPTH. Therefore, we will only study below the energy denominators of 
LFPTH diagrams along with the l/fc^-f actors for internal lines. 

We begin our analysis with Case B from Table [T] in Appendix \K\ x > A, illustrated in the 
diagram (B) in Fig. ^ Concentrating on the light-front energy denominators and 1//;; "'"-factors 
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(A) 




k+q q+r 




(B) 




Figure 9: The diagrams in LFPTH contributing to the SIDIS aniphtude in Fig. 3 
labels (A) and (B) denote the diagrams corresponding to Cases A and B in Table 1 



(A). Here 



for internal lines, we write with the help of Eq. (21) 
1 



(A;+)2 {p+ - k+) {k+ - r+) 



X 











1 






p- 


k\ 

k+ 


fci+A2 , . 




p~ + q- 


k+ 


-'"^+'' +^e\ 
p+-k+ ^ *^ 









1 




p- 


+ q- 




(fc±-rx)2 
k+-r+ 


fc2+A2 , . 1 



(69) 



Employing Eq. ( [2l| again, which implies the light-cone energy conservation condition for the 
diagrams in Fig. [91 



_ , „ (gl + rx)2 ri + \^ 

p +q = T — + -; — T 



(70) 



and using the scale hierarchy (34) to simplify the second denominator, we recast Eq. (69) as 

A(l-2;)(1- A) 



+ n2 



P^Q 



P, +62 



x - A + 



2gx-(a:^r_^-Afcx) 



+ ie 



((1 - A) A;x - (1 - x) rl)2 + A^ (x - A)^ 



with (cf. Eq. (30)) 



xy -xil-x)M^ >Q 



(71) 



(72) 



to impose proton stability (cf. (A7)). Since the diagram numerators in the LFPTH and in 



the above Feynman diagram case are the same, the argument from Sec. |2.2| about the need 
for a complex phase to generate the spin asymmetry still applies. Therefore we need to take 



an imaginary part of Eq. (71), which arises only from the second denominator, thus putting 
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the intermediate state involving the k + q and p — k hnes in Fig. M (B) on energy shelL The 



imaginary part of Eq. (71) is 



-ttA 



p+Q' 



P, +a2 



ik± - f±y 



X 



1 + 



2gl • (fc± -f±] 



A 



(73) 



where we also expanded the delta-function prefactor in powers of l/Q"^, which, among other 
things, put b = a. 



Comparing Eq. (73) to Eq. (28) (which was also the result of the diagram evaluation in 
Case B) or to Eq. (29), we see that the denominators of LFPTH give the same structure as the 
Feynman diagram calculation. 

To study Case A, x < A, we analyze the graphs in Fig. M (A). We easily observe that the 
second (right-panel) diagram in Fig. M (A) has no imaginary part and can thus be neglected 
for our purposes. The difference between the first (left-panel) diagram in Fig. M (A) and the 
graph in Fig.|9](B) is in the third energy denominator (corresponding to the latest intermediate 
state), which, together with the l/{r~^ — k^) factor, gives 



1 



P" + g" 



k+ 

-xA 



(rx-fc±)^ 
r+-k+ 



f^+A2 



+ ie 



(A - x) gl + A A?j 



xr_|_ 



{kl. - r±) 



(74) 



where, in the last step, we have used the delta-function from Eq. (73) common to the imaginary 



part of both diagrams along with the scale hierarchy (34). We see that the third intermediate 



state gives the same contribution to the imaginary parts of the diagrams in Fig. ^ (A) and (B) 



hence Eq. (73) is also valid in Case A. 



This completes our demonstration of the equivalence of the LFPTH calculation for the 
single-spin asymmetry in SIDIS to the Feynman diagram calculation. 



4.2 DY 

The LFPTH analysis of the Drell-Yan process proceeds along the lines similar to the SIDIS 
case. The LFPTH diagrams contributing to the phase-generating amplitude in Fig. IS] (A) are 



shown in Fig. 10 and are labeled (A), (B) and (C) according to the three non-trivial cases listed 
in Table [2] of Appendix |B} 

Starting with Case A from Table [2| < x < A, we write the contribution of the light-cone 
energy denominators and l/A;+-type factors for the diagram in Fig. 10 (A) as 



k+ (r+ - A;+) (g+ - k+) {p+ - A;^ 

1 



X 







P k+ r+-k+ p+-r+ ^ '*" 



p 


1 (9±-rx)2 
^ q+-r+ 


k+ 


q+-k+ 


rl+A2 , . 
p+-r+ 



(75) 



24 



M) 



q-r q-k q 



P P-k 




p-r 



(B> 




^\/\f\J\/\j 



q-r q-k 




(C) 




Figure 10: The diagrams in LFPTH contributing to the Drell-Yan aniphtude in Fig.p^(A). Here 
labels (A), (B) and (C) denote the diagrams corresponding to Cases A, B and C in Table [2] of 
Appendix iBl The labels on all the lines indicate the momentum flowing to the right. 



Using the kinematics in Eq. (47) leading to the light-cone energy conservation condition 

(g± - n_)2 



V 



q-i- — r 



q' + 



rj + X' 



we rewrite the imaginary part of Eq. (75|) in the following form 

-ttA 



P+Q^ 



B +a? 



S 



{ki. - r-l)2 



X 



in agreement with the factors in Eq. (53) and/or in Eq. (54). 



(76) 



1 + ^^^^^^lA 



(77) 



The analysis in Case B, A < x < /3, is slightly more involved. The denominators of all four 
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graphs in Fig. 10 (B) combine to give 



p- 



k — {p — k) +ie {q — r) — {q — k) — {k — r) +ie 



1 



p^ + {q — r)~ — {q — k) — {k, 
1 



{p-ky 



le 



-k- 



{q-ky 



le 



(P 



(k 



{p-ky 



le 



\p — k — {p 



k) +2e][(g — r) — {q — k) — {k — r) +ie] 
1 1 



te 



(7J 



{p — r) — [k — r) —{p — k)^+ie q — {q — k) —k 

The first term in the curly brackets of Eq. ( [78| , multiplied by the other two denominators, 
contains no imaginary part, since each energy denominator in that term represents a forbidden 
1 — )■ 2 decay or a forbidden 2 — )■ 1 merger. Only the second term in the curly brackets of 
Eq. (78), corresponding to the diagrams on the left of Fig. 10 (B), has an imaginary part. A 



simple calculation shows that this imaginary part (after it is multiplied by the l/fc^-type terms) 



is equal to Eq. (77), thus extending its validity into the A < x < /3 region in agreement with 



our Feynman diagram calculations. 

Finally noticing that the diagram in Fig. 10 (C) has no imaginary part, just like in the Feyn- 
man diagram case, we complete the demonstration of the equivalence of the LFPTH calculation 
for the single-spin asymmetry in DY to the Feynman diagram calculation. 



5 Conclusions 

In this paper we have calculated the single transverse spin asymmetries for the •y* +p^ —> q + X 
and g -|- p^ — ;• 7* + X processes in the model where the proton consists of a quark and a scalar 
diquark. We have shown explicitly that the SSAs arise from different cuts in the two processes 
in the Feynman diagram language (see Figs. Ill and |6]), corresponding to putting different energy 
denominators on energy shell in the LFPTH formalism. In spite of this difference, in the end of 
the calculation we get a simple sign-flip relation between spin asymmetries in the two processes. 



Eq. (62), in agreement with the arguments based on time-reversal anti-symmetry of the SSA 
jS]. The detailed calculation is consistent with the underlying dynamics of the lensing effect 
shown in Fig. [2j for QED. The flnal-state interaction in SIDIS is attractive, whereas the initial 
state interaction is repulsive in DY. Note that the Sivers effect is leading twist in Q^ even 
though the virtuality of the exchanged gluon which appears in the lensing effect is small. This 
is consistent with the OPE which is valid for small values of the ratios M'^/Q^ and r'j_/Q'^. 

The origin of the sign reversal at the diagrammatic level is discussed at the end of Sec. |3.2 
following Fig. [7j It appears from this discussion that the sign-flip between the SSAs in SIDIS 
and DY may only hold for the large-s and large-Q^ kinematics considered here and given by 



Eq. (19). While this was not checked explicitly in this work, it appears that the sign-flip relation 



(62) may not hold outside of this approximation, and may thus be destroyed by corrections to 



the transverse-momentum distribution (TMD) factorization used in \8\. More work is needed 
to investigate this further. 
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To summarize, we have confirmed the sign-flip relation (62) by an explicit diagrammatic 



calculation in a simple and robust model, in the process finding subtleties in the diagrammatic 



representation of (62) which were not known before. 



Acknowledgments 

Yu.K. is grateful to Leonard Gamberg, Jianwei Qiu, and Feng Yuan for discussions which 
inspired him to revisit the problem. The research of S.J.B. was supported in part by the 
Department of Energy contract DE-AC02-76SF00515; the research of D.S.H. is supported in 
part by the Korea Foundation for International Cooperation of Science & Technology (KICOS) 
and the Basic Science Research Programme through the National Research Foundation of Korea 
(2012-0002959); the research of Yu.K. and M.S. is sponsored in part by the U.S. Department 
of Energy under Grant No. DE-SC0004286; the research of I. S. is supported by Project Basal 
under Contract No. FB0821, and by the Fondecyt project 1100287. 

A Integration in the SIDIS Case 



Our goal in this Appendix is to evaluate the expression (27). There are four poles to the dk" 
integral, labeled below as (D - ®. Depending on the hierarchy of the longitudinal momentum 
fractions x and A, these poles may be located either above or below the real k~ axis. Since the 
outgoing quark and scalar are on-shell, we have (g+r)"'' = Ap'^ > and (p— r)+ = {1—A)p'^ > 
so that < A < 1. This allows us to write four distinct kinematic regimes in which to classify 
the poles: (x < < A < 1) , (0 < a; < A < 1), (0 < A < a; < 1), and (0 < A < 1 < x). 
The classification of the four pole locations as above or below the real k~ axis for each of these 
regimes is listed in Table [TJ 



(D 



Pole 



k- 
k- 
k- 
k- 



xp~' 



-q- + 



r + 
p- - 



{kj 



xp+ 



(a;-A)p+ 
{l-x)p+ 



Contribution: 



X < 



above 
above 
above 
above 







< X < A < 1 



below 
below 
above 
above 



Case A 



< A < X < 1 



below 


below 


below 


below 


below 


below 


above 


below 



Case B 



X > 1 







Table 1: Table classifying the pole locations of (27) as lying either above or below the Re/c 
axis. 



For X < or X > 1, all the poles fall on the same side of the Re k~ axis, so that we can close 
the contour in the other direction and get zero contribution. The physical region corresponds 
to < X < 1, and there are two distinct time-orderings of the diagram, x < A and x > A. We 
examine these two cases below. 



For Case A:0<x<A<l, we can close the contour in the upper half-plane, enclosing 
the poles (D and ®, as in Fig. 11 Let us consider the possible contributions to (27) from the 
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Imfc 




Refc" 




Figure 11: Illustration of the poles (left) and corresponding time-ordered diagram (right) of 



(27) for the kinematic regime Case A:0<x<A<l. We choose to close the contour in the 
upper half-plane, enclosing the poles ® and ®. Note that the placement of the poles is only 
schematic, indicating the sign of their imaginary part; the placement on the real axis has no 
significance. 



residue and imaginary parts of the various poles. 
• Resfd)] Im[(D]: Kinematically Prohibited 
This term would yield a contribution of 



-2tt Im 
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(Al) 



but the argument of the delta function is positive definite, so it cannot be satisfied; this cut is 
kinematically prohibited because it corresponds to a 2 —t- 1 massless, on-shell process. 

• Res[(i)] Im[®]: Proton Decay 



Similarly, this cut would yield a contribution of 
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(l-A)(A-a;)A;^ 
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All of the terms inside the delta-function are negative definite except for the first one, so we 
can impose the stability of the proton by requiring that 



A2-(1-A)M2>0. 
• Res[(D] Im[(D] + Res[®] Im[(D]: Kinematically Prohibited (Cancels) 
The term corresponding to Res[(D] Im[®] is 



(A3) 
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The sign of the ie argument of the imaginary part is ambiguous; this is typically a signature of 
a false pole. Whatever the sign of this term, it is exactly canceled by the Res[®] Im[®] term: 
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Thus this cut, which would correspond to a massless, on-shell 1 — )■ 2 process, is kinematically 
prohibited. 

• Res[®] Im[®]: Proton Decay 



This cut corresponds to proton decay through a different channel, yielding a contribution 



of 



X 



X (X S 





+271^ 
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(A6) 



-x{\^-{l- x)M^) - k] 



To prevent proton decay through this channel, we need to impose the slightly different condition 

A^ - (1 - x)M'^ > (A7) 
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Res[®] Im[(D]: Legal Cut 



This combination is the only legal cut of the four denominators that can be put on-shell 
simultaneously. This contribution is 
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+2ti^5 


p +q 
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{l—x)p+ xp+ 
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_ *^^ 1 

xp+ 




p~ 


' {l-x)p+ ' (A-x)p+ 



(A8) 



Expanding the argument of the delta function and keeping terms of order O {-q\ gives 
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kl + A^ (gl + kj 



(1 — x)p+ 



xp^ 
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1 + 2 



q± ■ (k± - f±] 



A 



(A9) 



The delta function sets s ^ A to leading order, but the singularity of the delta function only 
falls within the kinematic region of Case A, 0<a;<A<lif 

q± ■ {k± -f±) < 0, 

which is restricted to only half of the total phase space of the d'^k integral. As we will see, Case 
B complements this integral with the other half of the phase space. With this caveat, we can 
write a final expression for the imaginary part as 
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27r2AV 
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For Case B: < A < x < 1, we can close the contour in the upper half-plane, enclosing 



only the pole ®, as in Fig. 12 Again, we will consider the various contributions to (27) from 



the residue and imaginary part of the various poles. 



Im Ar 




Re fc- 




Figure 12: Illustration of the poles (left) and corresponding time-ordered diagram (right) of 



(27) for the kinematic regime Case B:0<A<x<l. We choose to close the contour in 



the upper half-plane, enclosing only the pole ®. Note that the placement of the poles is only 
schematic, indicating the sign of their imaginary part; the placement on the real axis has no 
significance. 
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Res[®] Im[(D]: Proton Decay 



This cut is exactly the same as the corresponding cut (A6) in Case A; it is unaffected by 



changing the sign of (a; — A). Thus the condition to prohibit proton decay through this channel 
is the same: A^ - (1 - x)M^ > 0. 

• Res[®] Im[(3)]; Kinematically Prohibited 

In this regime, it is explicitly impossible to perform this cut: 
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Since the argument is negative definite, this cut is kinematically prohibited. 
• Res[®] Im[(D]: Legal Cut 
Again, this is the only combination of propagators that can be put on shell simultaneously. 



The expression is the same as in (|A8|) but with x > A. This means that the delta function 



X 



r±) 



Q' 



A 



has its singularity within the kinematic window of Case B, < A < x < 1, if 

q± ■ {kj_ - fj_) > 0. 



Thus Case B gives rise to the same expression (AlO) at leading order, but with validity in 
the complementary region of the d'^k phase space; the final expression (AlO) is thus valid for 
all /c I . 



The expression (AlO) is illustrated in Fig. 4 in terms of the extra cut corresponding to 



putting the quark and the diquark propagators on mass shell. 



B Integration in the DY Case 



Here we want to evaluate the expression in ( 52 ) . In Table 2 we classify the four poles ® - ® of this 



expression as lying either above or below the Re k axis for the five distinct kinematic regimes: 
(x < < A < /3 < 1), (0 < a; < A < /3 < 1), (0 < A < a; < /3 < 1), (0 < A < /3 < s < 1), and 
(0 < A < /3 < 1 < x). As before all the poles lie to one side of the real axis unless < a; < 1, 
so there are three distinct cases to evaluate, each of which corresponds to a particular time- 
ordering of the diagram. We consider each of these cases below. 

For Case A:0<x<A</3<l, we choose to close the contour in the lower half-plane. 



enclosing only the pole ®, as shown in Fig. 13 Let us consider the possible contributions to 



(52) from the residue and imaginary parts of the various poles. 
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Table 2: Table classifying the pole locations of (52) as lying either above or below the Re/c 
axis. 
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Figure 13: Illustration of the poles (left) and corresponding time-ordered diagram (right) of 



(52) for the kinematic regime Case A:0<x<A</3<l. We choose to close the contour in 



the lower half-plane, enclosing only the pole ®. Note that the placement of the poles is only 
schematic, indicating the sign of their imaginary part; the placement on the real axis has no 
significance. 

• Res[®] Im[®]: Proton Decay 



This contribution would give 
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so we can prohibit proton decay through this channel by requiring that A^ 
• Res[®] Im[@]: Proton Decay 



;1-x)M2>0. 
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This contribution would give 
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xp+ 
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- 1 (fcx-rx)2l 
^ {A-x)p+ _ 




\kl 

xp+ 


y ^ (i3~x)p+ _ 




"1 
a;p+ ^ 


fci+A2 " 
(l-x)p+ 



(B2) 



oc 5 



X 



(A - x) (A^ - (1 - A)M2) + (1 - A)(A - x)kl + x(A - x)r=l 



+ x{l - A){k^ - r^] 



All of the momenta are positive definite, so we can prohibit proton decay through this channel 
by requiring that A^ - (1 - A)M^ > 0. 



Res[®] Im[@]: Legal Cut 



This corresponds to the only legal cut of the diagram as shown in Fig. 13; it is permitted 



because it corresponds to a 2 — )■ 1 process in which the two massless quarks become a single 
"massive" time-like photon with "mass" Q. Equivalently, we can recognize that the subsequent 
leptonic decay of the time-like virtual photon makes this cut correspond to a massless, on-shell 
2 — )■ 2 scattering process, which is allowed. This cut makes a contribution of 



X 





+27r'^6 
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xp+ 


q 


_ (fex-g±)2l 

"^ (f}-x)p+ _ 




\kl 

xp+ 




\—x)p+ 




"1 

[xp+ P^ + 


kl+x^ ' 

(l-x)p+_ 



(B3) 



and 



(k, 



q±} 



xp^ 



{/3 — x)p+ 



A(/3 - A)p^ 

0" 



X- 1 



— * 

q± ■ {kj 



r±; 



Q' 



A 



As usual, the 5-function sets a; ~ A, but the singularity only falls within the kinematic window 
of Case A (x < A) for q± ■ {k± — f±) < 0. As with DIS, this half of the (Pk phase space will be 
complemented by an equal contribution for Case B A < x < /3. Thus, the legal cut gives 



I 



2n^A{/3 - A)p+ 




s 


X - f 1 + 2 
V 


gx-(fc±-rx)A A 




Q' 


iki 
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- r~ 


1 (fcx-r±)2l 
"•" (A-a;)p+ 




"1 

[xp^ P~ + 


kl+\^ ' 

il-x)p+\ 



(B4) 



For Case B: < A < x < /3 < 1, we close the contour in the lower half-plane, enclosing the 
poles (D and (D, as shown in Fig. 14 Let us consider the possible contributions to (52) from 



the residue and imaginary parts of the various poles. 
• Res[®] Im[®]: Proton Decay 



The evaluation of this cut proceeds along exactly the same lines as in (Bl) of Case A; the 
proton stability condition is unaffected by changing the sign of (A — x). 
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Iinfc 




Refc- 




Figure 14: Illustration of the poles (left) and corresponding time-ordered diagram (right) of 



(52) for the kinematic regime Case B:0<A<x</3<l. We choose to close the contour 



in the lower half-plane, enclosing the poles ® and ®. Note that the placement of the poles is 
only schematic, indicating the sign of their imaginary part; the placement on the real axis has 
no significance. 

• Res[®] Im[(i)] + Res[(3)] Im[(D]: Proton Decay (Cancels) 



Evaluating Res[®] Im[(I)] would give a contribution of 



Xi 



-h27rlm 



(x — A)p+ 



xp+ y "•" (fs-x)p+ 



xp+ " (1— x)p+ 








±271^6 


xp+ 
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x-gx)^! 

3-x)p+ \ 
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fc2^+A2 " 
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(B5) 



where the sign ambiguity of the it components indicates the presence of a false pole. Whatever 



the sign of (B5), it is exactly canceled by the contribution of Res[(i)] Im[®]: 



+27rlm 



(fc^-rj^)^-ie 
(a; — A)p+ 



„- „- 1 (fcx-rx)2 1 (fcx-gx)2" 
/ y "^ (x-A)p+ 1 (/3-x)p+ _ 




[' P ^ (x-A)p+ 1 (l-x)p+J 



^2Tx'^b 



r -|- 



(fcx-rx)' 
(a::-A)p+ 



xp"f 



r- _ g- + (fcx-rx)2 ^ (fcx-gx)^ 



P "•" (x-A)p+ "^ (l-z)p+ 



(B6) 



(x-A)p+ ' (/3-x)p+ 

Thus Xi -I-X2 = 0, so that proton decay by this channel is automatically prohibited for Case B. 
• Res [(3)] Im[(D]: Kinematically Prohibited 
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This contribution would be 



X 



+2ti^5 


r — q 


"•" {x-A)p+ ' (P-x)p+ 




Y ^ ix-A)p+ 


xp+ 




^ P ^ {x-A)p+ + 
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(B7) 



oc 6 



X 



A)u + i(3-x)r^-i(3-A)kj 



since the argument is positive definite, this cut is kinematically prohibited, as it corresponds 
to a 1 — )■ 2 niassless process. 

• Res[@] Ini[®]: Kinematically Prohibited 



This contribution would be 



I 



+27lH 


r — 


p 


- {fcx-rx)2 fci+A2 1 
^ {x-A)p+ ' {l-x)p+\ 




L- 1 {k±-r±y 
' ^ {x-A)p+ 


k] 1 
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(B8) 



oc 6 



X - A)^A" + (1 - x)f^ - (1 - A)A;j 



since the argument is positive definite, this cut is kinematically prohibited. 
• Res[®] Im[(D]: Legal Cut 



Again, this is the only legal cut of the diagram in Fig. 14 The expression is the same as in 



(|B3|) from Case A, but with x > A. This means that the delta function 

q± ■ {k± - r±] 



X- 1 + 2 



Q' 



A 



has its singularity within the kinematic window of Case B,0<A<a:</3<l, if/c^- (qj_ —r±)> 
0. Hence we again recover (B4), but with validity in the other half of the (Pk phase space. 



Cases A and B thus complement each other, and we will show that Case C does not make any 



contribution to (52). 



For Case C: < A < /5 < z < 1, we choose to close the contour in the upper half plane. 



enclosing only the pole ®, as illustrated in Fig. 15 We demonstrate below that there is no 
viable cut for this time-ordering of the process. 

• Res[®] Im[®]: Proton Decay 



This cut is the same as in (Bl) of Case A; closing the contour in the other direction does not 



affect the overall sign of the contribution once the imaginary part is taken, and changing the 
signs oi P — X and A — x does not affect the proton stability condition. 

• Res[®] Im[@]: Kinematically Prohibited 
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Im/c 




ReA;- 




Figure 15: Illustration of the poles (left) and corresponding time-ordered diagram (right) of 



(52) for the kinematic regime Case C:0<A</3<a;<l. We choose to close the contour in 
the upper half-plane, enclosing only the pole ®. Note that the placement of the poles is only 
schematic, indicating the sign of their imaginary part; the placement on the real axis has no 
significance. 



This cut is the same as in (B8) of Case B; closing the contour in the other direction does not 



affect the overall sign, and changing the sign of /3 — x does not affect the argument of the delta 
function. Hence this process is also kinematically forbidden. 

• Res[®] Im[(D]: Kinematically Prohibited 



This is the only new cut that requires explicit calculation. This contribution would be 



+2tx^5 



X 



P -Q 



{l-x)p+ {x-l3)p+ 



P - 



(l-x)p+ 






p — r 



fcj+AZ (fcx-rl)2 



{l-x)p+ {x-A)p+ 



(B9) 



oc 6 



{(3 - A)(x - (3){x - A)\^ - (1 - x){x - /3)(1 - /?) r^ - 



1- A 



q± 



- (l-A)(/3-A)(l-/3)(^fc^-i-|g-l 
Since the argument is negative definite, this process is kinematically forbidden; this cut would 



not only correspond to proton decay in the lower half of the diagram in Fig. 15, but also a 
kinematically prohibited 3 — )■ 1 process in the upper half. Thus we have shown that there is 
no viable cut of the diagram for the kinematics of Case C, and this case makes no contribution 
to the asymmetry. Therefore ( |B4[ ) gives the complete expression for the imaginary part and is 
our final result. 



(B4) is illustrated in Fig. Mm the text by the second (shorter) cut putting the quark and 



anti-quark propagators on mass shell. 
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